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Abstract 

We estimate the quantum state of a light beam from results of quantum homodyne measurements per- 
formed on identically prepared pulses. The state is represented through the Wigner function, a "quasi- 
probability density" on which may take negative values and must respect intrinsic positivity constraints 
imposed by quantum physics. The data consists of n i.i.d. observations from a probability density equal to 
the Radon transform of the Wigner function. We construct an estimator for the Wigner function, and prove 
that it is minimax efficient for the pointwise risk over a class of infinitely differentiable functions. A similar 
result was previously derived by Cavalier in the context of positron emission tomography. Our work extends 
this result to the space of smooth Wigner functions, which is the relevant parameter space for quantum 
homodyne tomography. 

1 Introduction 

The phenomena occurring at the interface between the microscopic and the macroscopic worlds have an intrinsic 
probabihstic nature. When measuring properties of atoms and laser pulses we obtain a random result whose 
probability distribution is determined by the state, or preparation of the quantum system. For example, if 
we count the number of photons coming from a laser source we observe a Poisson distributed random variable 
with mean equal to the intensity of the laser. The statistical inverse problem of inferring the state from results 
of measurements on many identically prepared systems, is called quantum state estimation. Recently it has 
become possible to apply such a method to the reconstruction of the quantum state of a light beam. The 
measurement technique is called Quantum Homodyne Tomography [22] and is used to confirm the creation of 
new and exotic quantum states of light such as squeezed states [3] , single-photon-added coherent states [H] and 
Schrodinger cat states [20] . As experiments become more and more complicated, the costs - in terms of money 
and time - of running a measurement rise, and one needs to apply more sophisticated statistical techniques to 
reconstruct the state from a limited number of data. This paper makes a step in this direction by providing 
minimax convergence rates for a class of physical states. 

The object to be estimated is a real function of two variables called the Wigner function [23J, which can be seen 
as a joint density of the electric and magnetic fields of the laser beam. However, since in quantum mechanics 
we cannot measure both electric and magnetic fields simultaneously, this function is in general not a probability 
density but has many features in common with the latter, for example the marginals along any direction are 
bona-fide probability densities. In quantum optics the Wigner function is a preferred representation of the 
quantum state [5] because many interesting quantum patterns such as squeezing or oscillation between negative 
and positive values, can be easily identified from the shape of the function. The estimation methods used by 
the physicists involve a number of ad-hoc approximations, binnings and truncations, making it difhcult to verify 
the reliability of the procedure. Moreover, the quantum features in which the experimenter is interested may 
be washed out in the resulting estimator. 
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From the statistical viewpoint, we deal with an iU posed inverse problem which is closely related to the problem 
of estimating a bivariate probability density in Positron Emission Tomography. In both cases the parameter 
of interest is a density and the data consists of independent identically distributed samples from the Radon 
Transform 7J of that density, with uniformly distributed angles. However as we mentioned above, the Wigner 
function is not necessarily positive but must satisfy other positivity conditions dictated by the laws of quantum 
mechanics. This means that some of the existing statistical results are not directly applicable in this case, and 
the problem of estimating the Wigner function should be studied separately in order to identify the specific 
"quantum features" that could be exploited in designing new estimators. The problem of state estimation in 
Quantum Homodyne Tomography was investigated in the non-parametric setting in [I]. The paper provides 
sufficient conditions for consistency of various estimators for two different representations of the state, namely 
the density matrix and the Wigner function. The problem of estimating the Wigner function was further 
investigated in [4j in a set-up which takes into account the detection losses occurring in the measurement, 
leading to an additional additive Gaussian noise. Although both [4J and the present paper consider a family 
of very smooth Wigner functions, the estimation techniques are very different. In |4] the bias dominates the 
variance due to the presence of the Gaussian noise, and for deriving the lower bound it is enough to consider a 
"worst family" consisting of just two states. In this paper, the variance is dominating and for the lower bound 
we need to consider a one parameter family of states. The techniques that we use here are very similar to the 
one of Cavalier [5] who considers the same tomography problem in the context of smooth probability densities 
rather than Wigner functions. For the problem of Positron Emission Tomography and the related regression 
problem of X-ray Tomography we refer to [161 1171 E] the references therein. For an introduction to Quantum 
Statistical Inference we refer to [H [T] . 

In Section [2] we give a short introduction to quantum mechanics, present some properties of Wigner functions 
and the relation to the density matrix which will play an important role in the proof of the lower bound. The 
statistical set-up and the main results are presented in Section [31 Following [^ , we identify a class of very 
smooth Wigner functions 



where p and L are positive constants. From the physical point of view, all the states which have been produced 
in the lab up to date belong to such a class, and a more detailed argument as to why this assumption is realistic, 
can be found in [4^. 

We consider a family of estimators depending on a bandwidth which is chosen according to the parameters of 
the class. The upper bound for the pointwise risk is proven in Theorem [T] and has the same almost parametric 
expression as the bound derived in [5]. For the lower bound we consider a family of Wigner functions which 
is different from the worst parametric family of probability densities of [S]. The latter cannot be used in our 
situation since it does not correspond to Wigner functions of physical states. Thus the main novelty of the 
paper is the derivation of the lower bound in the physical context of Quantum Homodyne Tomography rather 
than that of Positron Emission Tomography. 

The technical Lemmas are grouped in Section [Jj 

2 Quantum Homodyne Tomography and the Wigner function 

In this section we briefiy present some basic notions of quantum mechanics, the mathematical set-up of Quantum 
Homodyne Tomography, and some properties of Wigner functions. More information on the general set-up of 
quantum statistical inference can be found in the review paper |2j and the textbooks |13j and |14| . Quantum 
Homodyne Tomography is discussed in detail in [U [3| . 

The state of a quantum system encodes all necessary information for computing the probability distribution 
of results of any given measurement. Mathematically, the state is described by a density matrix, which is a 
compact operator p on a complex Hilbert space Ti. having the following properties: 

1. Selfadjoint: p = p* , where p* is the adjoint of p. 




(1.1) 
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2. Positive: p > 0, or equivalently {ip, pip) > for all tp gH. 

3. Trace one: Tr(p) = 1. 

Positivity implies that the eigenvalues of p are all nonegative and by the last property, they sum up to one. 
Notice that the above requirements parallel those of defining probability densities. 

We have the following diagonal form 

dim?^ 

P=Y1 ^^P' (2.1) 

i=l 

where pi is the projection onto the one dimensional space generated by the eigenvector G 7i of p corresponding 
to the eigenvalue Aj, i.e., pci = XiCi. With respect to a fixed orthonormal basis {4'i}i>i in H, the operator p 
can be represented as a matrix with elements pij = (ipijpipj). 

Let us consider now the following problem. We are given a quantum system prepared in an unknown state 
p and wc would like to determine p. In order to obtain information about the system we have to measure 
its properties. The laws of quantum mechanics say that for any given measurement with space of outcomes 
given by the measure space {X, S), the result of the measurement performed on a system prepared in state p is 
random and has probability distribution Fp over {Cl, S) such that the map 

p ^ ¥p, (2.2) 

is affine, i.e. it maps convex combinations of states into the corresponding convex combination of probability 
distributions. This has a natural interpretation: a system can be prepared in a mixture Xpi + (1 — A)p2 of states 
by randomly choosing the preparation procedure according to the individual state pi with probability A and p2 
with probability 1 — A. The distribution of the results will then refiect this randomized preparation as well. 

The most common measurement is that of an observable such as energy, position, spin, etc. Any given observ- 

ables is described by some selfadjoint operator X on the Hilbert space H and we suppose for simplicity that it 
has a discrete spectrum, that is, it can be written in the diagonal form 

dim?i 

X = ^ x„P„. (2.3) 

i=l 

with Xa G K and Pq one dimensional projections onto the eigenvectors of X. The result X of the measurement 
of the observable X for a preparation given by the state p, is a random element of the set Q = {xi,X2, . . . } of 
eigenvalues of X and has the probability distribution 

¥p[X = Xa]=TT{Pap). (2.4) 

This measurement will give statistical information about the diagonal elements of the density matrix p with 
respect to the eigenbasis of X, and it suggests that in order to estimate all matrix elements of p one would 
have to probe the system from a number of "directions" by performing different measurements on identically 
prepared systems. 

This can be generalized to the case of infinite dimensional Hilbert space, and measurements with outcomes in 
arbitrary measure spaces. In the next section we will see that an infinite dimensional density matrix can be 
estimated by measuring a randomly chosen observable from a continuous family of non-commuting observables. 

2.1 Quantum homodyne tomography and the Wigner function 

An important example of quantum system is the monochromatic light in a cavity, described by density matrices 
on the Hilbert space of complex valued square integrable functions on the real line, L2{M.). A distinguished 
orthonormal basis of this space is given by the vectors 

Vfe(ar) = ^^fe(ar)e-^'/^ fc = 0,1,2..., (2.5) 
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where Hk are the Hermite polynomials normalized such that tpk (x) is a unit vector representing the pure state 
of k photons inside the cavity. We will denote the matrix elements of p with respect to this basis by pij. Notice 
that the diagonal of the density matrix is a probability distribution over the nonnegative numbers pk — Pk,k- 
This is the probability distribution of results when measuring the number of photons in the cavity prepared in 
state p. Clearly this distribution does not contain information about the off-diagonal elements of p thus it is 
not sufficient for identifying the state of the system. This is a typical situation in state estimation and one has 
to devise experiments in which the systems are looked at subsequently from "different directions" , a broadly 
described methodology which in the physics literature goes by the name of quantum tomography. Quantum 
Homodyne Tomography is one such measurement method which is frequently used in quantum optics at the 
estimation of the quantum state of light [31 [201 121] ■ We skip the measurement set-up which is described in detail 
in [II S] and present the statistical problem associated to this measurement. 

We observe {Xi, $1), . . . , (X„, $„), i.i.d. random variables with values in R x [0, tt] and distribution Pp whose 
density with respect to the measure Tr~^ d(j) x dx is given by 

00 

= J2 P,,kMxHj{x)e-'^'-''^^. (2.6) 

3,k=0 

Since p is a positive definite matrix of trace 1, and ipj form an orthonormal basis, it follows that pp is a 
probability density: real, nonnegative, integrates to 1. The data (Xi, ^i), £ — 1, . . . ,n, come from independent 
measurements on identically prepared pulses of light escaping from the cavity, whose state is completely encoded 
in the matrix p. 

For each of the systems independently, we repeat the following experimental procedure: we first choose the 
angle $ uniformly distributed over [0, tt] and then measure a certain observable called quadrature, obtaining 
a real valued result with probability density Pp{x,(j)). The quadrature is defined as the linear combination 
X^ := cos(/)Q -I- sin0P, where Q and P are the electric and magnetic fields of the light beam given by the 
selfadjoint operators 

Q^(a;) = xil}{x), and PV'(a;) = —i-f-. 

dx 

The characteristic functions of these densities can be put together to define a function of two variables 

Wp{u,v) := Tr(pexp(-iiX^)) = Ti[pp{-A)]{t), (2.7) 

where we have used the polar coordinates (u, w) = {t cos (j), t sin (p) , and Ti is the Fourier transform with respect 
to the first variable, for fixed </>. Note that our convention for defining the Fourier transform and its inverse are 
the following 

/oo -1 /"OO 

/(x)e-"*dx, T-\g]{x) = - / g(i)e"*di. 
-00 ^'^ J —00 

Equivalently, we can write 

Wp{u,v) =Tr(pexp(-mQ- wP)), (2.8) 

which resembles a characteristic function of a bivariate probability density, namely the joint density of Q and 
P. However, since the operators Q and P^ do not commute with each other, we cannot speak of their joint 
probability distribution and the function Wp{u,v) is in general not the characteristic function of a probability 
density but rather of the so called Wigner function 

Wp{q,p):=T^\Wp]{q,p), (2.9) 

a "quasi-distribution" which may take negative values but whose marginals are bona-fide probability densities. 
As we will see below, the Wigner function Wp is in one to one correspondence with the density matrix p, and 
in quantum optics one frequently uses the Wigner function as an alternative representation of the quantum 
state, having the advantage that it illustrates important "quantum features" such as squeezing and negative 
oscillations. From (|2.7p and (|2.9p we deduce that the probability density of the data Pp{x,4>) is the Radon 
transform of the Wigner function 

/oo 
Wp{x COS (I) — t8in4>,xsm4> + tcos (j))dt. 
-00 
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adding Quantum Homodync Tomography to a long the list of applications ranging from computerized tomog- 
raphy to astronomy and geophysics jTj . 

Another important feature of the Wigner function is that it can be used as a computational tool: for any 
selfadjoint operator X there exists a function Wx from M.^ to M such that the expectation of X is given by 

Tr(Xp) = 2^// W^{q,p)Wp{q,p)dqdp. (2.10) 

In particular, the correspondence between the density matrix p and the Wigner function Wp is an L2 isometry 
up to a constant 

\\Wp-Wr\\l:^ ff \Wp{q,p)~Wr{q,p)\'dpdq^^ £ - r.,,f • (2.11) 



The space of Wigner function has an overlap with that of probability distributions. For example, all Gaussian 
densities which are bounded from above by I/tt are Wigner functions, while the rest do not correspond to physical 
states. This is due to the celebrated Heisenberg's uncertainty relations which say that the non-commuting 
observables P and Q cannot have probability distributions such that the product of their variances is smaller 
than J . In general, a Wigner function cannot be too "peaked" : 

\Wp{q,p)\<-, forall(q,p)eR2_ (2.12) 



Some examples of quantum states which can be created in laboratory are given in Table 1 of pj. Typically, 
the corresponding Wigner functions have a Gaussian tail but need not be positive. For example the state of 
one-photon in the cavity is described by the density matrix with pi.i = 1 and all other elements zero which is 
equal to the orthogonal projection onto the vector ^i. The corresponding Wigner function is 

W{q,p) = -{2q^ + 2p2 _ 1) exp(-g2 _ p^). 



In conclusion, although we deal with a problem which is similar to that of Positron Emission Tomography, the 
parameter space is different from the space of probability densities and special techniques have to be developed 
for this situation. 



3 The main results 

Our problem is that of estimating the Wigner function Wp{z), defined on the plane z = {q,p). In order to prove 
rates of convergence some restrictions are necessary to be imposed to the class. We consider the class yV(/3,L) 
of Wigner functions which are continuous and whose Fourier transform satisfy 

( \W{w)\^ eM2(3\w\)dw < L (3.1) 

for [3 and L positive constants. This condition implies that the function to be estimated is very smooth. Such 
classes appeared in the statistical literature in [15j , and we subsequently used in the context of density estimation 
[llj . functional estimation, regression problems [lOj . and tomography 5 . In 4J we have argued that from the 
physical point of view it is natural to assume that the Wigner function of a state which can be created in the 
lab belongs to such a class. 

We use a kernel-type estimator based on the following function called a band-limited filter 

= ^e-. (3.2) 
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This filter has already been used in the context of tomography [1^1 [13 [S] , and its Fourier transform is 



Ksjt) ^ -\t\IsAt), 



(3.3) 



where Is„ is the indicator function of {t : \t\ < 1/Sn}. The estimator we use is 

— 1 " 

W„(z) = -yKsJ[z,^,]-X,) 
n ^ — ' 

i=l 

with i.i.d. observations {Xi, for i = 1, . . . , n, with density Pp{x, (f>) — TZ[Wp\{x, (p). 
Following [IH] we define the dual operator 7?.^ on Li(M x [0, tt]) by 



Then 



n*[h]{z 



n*n[w]{z) 



h{[z,cj)],(j)) dcf). 



(3.4) 



(3.5) 



(3.6) 



represents the integrals of W over all lines passing through the point x. Note that in general R^Tl[W]{z) > 
for all Wigner functions W and all z e M^, and the number states V'fc with k odd have the property that 
TZ'^Tl[W](Q) = 0. In it is assumed that the probability distributions / to be estimated are strictly positive 
which implies that TZ'^TZ[f]{z) > for all z. For the upper bound we will assume that the latter condition 
holds. 



3.1 The upper Bound 

Theorem 1 For any W G VV(/3, L) and any fixed z G M-^ such that R'^TZ[W]{z) > we have as n ^ oo, 



E 



(Wniz) - Wiz)f] = C*TZ*n[W]{z) X + o(l)) 



(3.7) 



where C* 



3(47i-/3)3 • 



Proof. We will provide only the main steps of the proofs pointing out where the assumption on the class of 
Wigner function plays a role. For a more detailed proof of the bounds for the class of probability densities in 
A{P, L) we refer to [5]. The risk can be decomposed in two parts, the bias and the stochastic part 



E 



(W^„(z) - W{z)f\ = (E(M/„(z)) - W{z)) + E [{W,,{z) - E(M^„(z))) 
■.^hl{z) + al{z). 



(3.8) 



On one hand, using property (|3.3p of the kernel and the inverse Fourier transform, the bias can be written as 

U^) = J^j ^(«^)^ (l^l > ^) e"'^"'''^ dw. (3.9) 



By using Cauchy-Schwarz and the fact that W E W{f3, L) we get 

bl{z)<^e''^f''^{l + o{l)), 

On 

as 6n 0. With the choice l/(5„ ~ logn/(2/3) the bias upper bound becomes 



bl{z)<c^^{l + o{l)), 
n 



as n — > oo. 



(3.10) 



(3.11) 
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On the other hand, the variance is equal to 



aliz) ^-Yi^tKsA^,^] - X) = 
n 



E[k2 --(E[if,„([z,<f>]-X)])^ (3.12) 



rt n 



where (X,^) is a random variable with probability density pp{x,(f>) = TZ[Wp]{x,4>). The second term can be 
bounded as follows 

i (E[KsA[z,'^>] - X)]f < ttA^ I hM)e-^^^'"Uw X [ \Wiw)\^e^^^^^ dw = Oi-), (3.13) 

n [liryn J J n 

with 0(i) uniformly with respect to G W(/3, L). 
The first term is 

E [Kl{[z, $] - X)] = f I Kl{[z, 4>] - v)Pp{y, cf) dcPdy. (3.14) 

Now denote 

G(u) = 3 ( / rcos(ur)dr) , (3.15) 

and let Gsiu) = {l/5)G{u/5). We have 







3 . .1 



2 / 1 \ 3 
TT^ / 1 ^ 



Kl{[z,cj>]~y)pp{y,4>)dy^^^^-J (G^,, * 7^[VF](•, 0))([z, 0]). (3.16) 
Using [S] Lemma 4, we have that as (5 ^ 0, 

r 

{Gs * n[W]{-,<l>)) d(p = ■R*n[W]{z){\ + o(l)) + 0(5i/3). 
With this the second term of the variance can be written as 







E [KUz, ^ - X)] = ^ {£)n*nW]iz)il + oil)) + O ( (£ 
as 5n 0. Thus as n ^ cx3 



a^(z) < G*7^#7^[VF](z)^i^(l + o(l)) + O f ^ 

n \ n 

^ fG*7^#7^[M^](^)(l + o(l)) + o(iog-i 



□ 



If i?#7^[M^](z) > then we obtain the claimed constant. Notice that if R*n[W]iz) = then the convergence 
1 ^ 

rate is faster than °^ " . 



3.2 The lower bound 

In order to prove a lower bound result we consider the slightly modified class of Wigner functions 

W(/9, L, an) ^{We W(/3, L) : n*n[W]iz) > a„}, (3.17) 
for a sequence a„ such that lim„^oo ctn — and lim„^oo(a?i(logn)^/'^) = oo. Let us denote 

r„(W,z) = (^C*n*n[W]{z)^^^^^^ ' . (3.18) 
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Theorem 2 For a fixed z e M^, we have 



liminf inf sup ] 

n^oo WeW(/3,L,Q„) 



W^{z)-W{z) 
rn{W,z) 



> 1 



where infrr; denotes the infimum over all estimators ofW{z). 



Proof. The proof is based on the standard procedure of building a hardest parametric subfamily for the class 
W(/3, L) of Wigner functions of the form 

Wc^Wa^ + cga, (3.19) 
where c is a parameter in a neighborhood of the origin, Wa and ga are functions to be defined shortly. 

The essential point of the proof is that the family of probability densities /o + cga used in ,5,, is not always 
contained in our parameter space consisting of Wigner functions. For illustration we will show that for some 
parameters (3, the function /o defined in equation (29) of [S] is not the Wigner function of a quantum state. 
Indeed, suppose for the moment that this was the case, i.e. fo = Wp for some state p. Then by the rotation 
symmetry of /o, the density matrix p must be diagonal and 

oo 

Pp{x,(j)) = ^pfe,feVfe(a;). 

Using the inequality [S] ||'0a;||oo < k, where fc is a constant whose value is slightly bigger that 1, and the fact 
that J2 Pk,k = 1, we find that ||pp||oo < k. However, the Radon transform of the Wigner fmiction Wp is 

Pp{x,4>)^n[Wp]{x,4>)^n[h]{x,4>) = 

which violates the above bound for /3 < l/(7rfc). 



We thus define a parametric subfamily of 'W{(3,L) which is a suitable modification of the family considered in 
[5] in order to cope with this problem. 



Construction ofWa- Consider the Mehler formula, (see [8], 10.13.22) 



y z" ^\ g„(x)^e-"' = , ^ exp f-x^L^] . (3.20) 



n—Q 



Integrating both terms with fa{z) — a{l — z)" we get 

" 



n=0 

The Fourier transform of pa is 



p„(x,0) :=f;^„(x)2 f Uz)z^dz^ /"_^(f)_expf-a;2^ j (3 2i) 

n Jo 



l + Z 



W^{w) T[p^]{w) = exp (^-\w\^A±^^ dz (3.22) 

Notice that the normalization condition J pct = 1 is equivalent to Wq(0) = 1 which is satisfied for the chosen 
functions /„, thus Pa is a probability density corresponding to a diagonal density matrix with elements 

1 

k 



Plk= / ^'fcdz. (3.23) 
Jo 

We denote by Wa the Wigner function whose Fourier transform is defined in equation (|3.22p with a > a 
parameter to be fixed later. This function is considered in a more general form in [4], and corresponds to 
for e = . 
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Construction of ga- Let [5] 



Ha{s) 



1 



■ cos(sr) c?r, 



s G 



(3.24) 



(27r)2 Jo 1 + a-isinh2/3r 

where a > is a parameter which wiU depend on n as a = a„ = rt'' with < 77 < 1. The Fourier transform of 
ga is 



Let 



5a (z) 



47r 1 + a-1 sinh^ 

cos((2;, w)) dw, 



w 



2(27r)3 7r2 l + a-isinh2/3|u;| 



and its Fourier transform 



9aiw) = 



1 



Now, let us consider the family 
where the real parameter c satisfies 



47rl + a-isinh'^/3|u;| 

q 



w € 



\C\ < Ca 



(3.25) 
(3.26) 
(3.27) 

(3.28) 
(3.29) 



V^(loga)3/2' 

with q > sufficiently small. By translating with z in we obtain our hardest family for pointwise estimation 
at the point z 

W,^O^W,iC-z). (3.30) 

We will check now that belongs to the class W{P,L) for an appropriate choice of a in Wa, which means 
that is a Wigner function and 



1 

4^ 



W[:{w) 



47r2 



Wc{w) 



(3.31) 



By Lemma [5] we have € W(/3, L/A) for a small enough a > 0, and from Lemma 5 of |5j we have that 



(3.32) 



for all c < Ca. The last two conditions together imply (|3.3ip . 

Furthermore, has to be a Wigner function. As translations in the plane transform Wigner functions into 
Wigner functions, we need only to show this for Wc- This means that there exist a family of density matrices 
Pc such that their corresponding Wigner functions are Wc- The invariance of Wc under rotations in the plane 
translates into the fact that p'^ has all off-diagonal elements equal to zero, and thus we only need to show that 
all its diagonal elements are positive and add up to one. The relation between the diagonal matrix elements 
and the Wigner function is [TS] 



Pk,k 



1 

2^ 



e\^\"^^Lki\w\y2)Wc{w)dw, 



where Lk are the Laguerre polynomials. By linearity we have 

Pfc,fc — Pk.k + CTfe^fc 

where 



'k,k 



and 



te-'"'^Lk{t'/2)Ha{t)dt, 



pt,=a / z\l-zr. 



(3.33) 

(3.34) 
(3.35) 

(3.36) 
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Corroborating the result shown m Lemma |4] 



as a, fc — *■ CX3, with that of Lemma 2 in ^ 

we conclude that if a < 1/4, then p^, > for all fc > and |c| < Ca for a sufhciently large. 

Now we can use the fact that for the family of translated functions, as defined in eq. (|3.30p . TZ'^TZ[W^]{z) — 
n*7i[W]{0). Indeed, using we get 



n[W']{x,(j>) 



1 
1 

2^ 



{t cos (j), t sm(j))e'*'' dt 

W{t cos (/), t sin 0)e-**[^'*le'*^ dt 

— /vF(tcos(/),tsin</))e**(^-[^''^l)dt 
27r J 

n[w]{x-[z,cj,],^). 



(3.37) 



Now, from definition of 7?.* we get 



27r 



7^[VK](o,?!))^^0 = n*n[w]io) 



Thus, := 7^#7^[VK^^](z) = 7^#7^[T4^„](0) + c7^#7^[ga](0) for any z. Given in our case and 5a are 

invariant under rotations we obtain 



n*n[Wa]{o) 



Pq(0, 0) dcf) 



1 (i_^)a-l/2 



> 0. 



/o (1 + ^)1/2 

For the second term, using eq. p.24[) . Lemma 5 in [5], and definition of Cq 

sup |c7^#7^[5a](0)| = sup |27rci/a(0)| = o(l) as a^oo. 

\c\<Ca \c\<Ca 



(3.38) 



(3.39) 



We conclude that Rc — Ro{l + o(l)) > q;„, for n — > oo and thus, for a large enough, G yV(/3, L,a„). 
Moreover, from p.lSp . 



r„(Ty,^z)2=r„(Ty„, 0)2(1 + 0(1)) = C*i?o^i^^(l + 0(1)). 



(3.40) 



The rest of the proof is based on the Van Trees inequality and follows along the lines of f5]. The main 
difference is in the proof of Lemma [6] where the Fisher information of the family of densities ^'R-\W^] is 
approximated. Take a continuously differentiable probability density, Ao(c), defined on the interval [—1, 1], such 
that Ao(— 1) — Ao(l) — 0, with a finite Fisher information /q. The new density A(c) = Aa(c) = C^^\o{C^^c) is 
a prior density with finite Fisher information /(A) = IqC^^. Finally let us define /(c) the Fisher information 
of the family of densities ^TZ[W^]. Using the Van Trees inequality 



> inf 



inf sup 

W„ WeW{(3,L,a„) 

Ca 

Ew 



(Wn{z)~Wiz) 



> inf sup Ew^ 

W„ \c\<Ca 



(w„{z)~W^iz) 



(Wn{z)-W!iz) 



Aq (c)dc > 



iSX {dW,Hz)/dc)Xa{c)dcf 

?^/fc„^W^a(c) dc + /(A) 



(3.41) 
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Now, from equation (|3?26| and Lemma 5 of [5], we get that dW^{z)/dc = ga{0) = C*{logaf{l + o(l)). By 
inserting the expression of /(c) from Lemma |6l and using p.40p and a = a„ = we get the lower bound 



inf sup E 

W„ \c\<Ca 



Wn{z)-W,^{z) 



> 



in{dw,H^)/dc)\cdcy 



(r„(W^,^z)(l + o(l)))- 



nJ_^J{c)Xcdc + I{\) 

(C*773(iogn)3)2r„(W^„,0)-2(l + o(l)) 



C*ifRo ^n(logn)3(l + 0(1)) + q-'^Io-q'^ni{\ognf ' 
By letting n ^ 00 followed by 77 — > 1 we finally obtain 



liminf inf sup E 



Wn{z)-W{z) 
r{W, z) 



4 Technical Lemmas 



> 1. 



□ 



Lemma 1 Let W € VV(/3, L). Then the following inequalities hold 

\W{uj)\<l, VcjgM^ 
\W{w)\dw < Q, 

sup n*n[w]{z) <i + 



W<^W(I3.L) 

where Q is a constant depending only on (3 and L. 



27r' 



(4.1) 
(4.2) 

(4.3) 



Proof. As = Wp for some density matrix p we use the definition of Wp to obtain 

\Wp{u,v)\ = |Tr (pe-»Q-™P) | < ||p||i = 1, 

where we have used the normalization of the density matrix p and the fact that \\U\\ = 1 for any unitary operator 
U. The inequality (|4.2p is a direct consequence of the definition of W(/3, L) and the Cauchy-Schwarz inequality. 
Now using (12. 7p and the previous inequalities we get 



\TZ*TZ[W]iz)\ 



2Tr 



n[w]i[z,4>],^)d^ 



1 



27r 



— / d(/) j dtW{t cos (j),t sin (f>)e 
1 

— / d(p I dt\W{t cos (j),t sin < 
27'' Jo 



it[z,(p] 



< 



1 

2^ 
1 

2jr 



{Q + 2tt). 



duj\Wiuj)\ 



d(j) [ dt\W{t cos (j),t sin < 



□ 



Lemma 2 For all W G yV(/3, L) we have 

\n[W]ix,^)-n[W]{y,(f>)\ d(l)<^\x-y\ 
with Q the constant depending on (3 and L defined in Lemmas 
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Proof. By (|2.7p and Lemma [T] we have 

r 

\TZ[W]{x,cl>)--n[W]{y,cl))\ dcf,< 



where we have used that |e**^ — e'*^| < 2|t||a:: — ?/|. 



^ - 27r 



< - y\, 



□ 



Lemma 3 For all < a < 1 and \x\ > 1 there exist constants c, C depending on a, such that 

Moreover there exist constants Ci and C2 such that the first two derivatives of pa satisfy \p'aix)\ < Cix^^^^^"\ 
and <C2X-(3+2"). 



Proof We have 



Pa{x,(l)) 



(1 + ^)1/2 



exp —x 



which by the change of variables u ^ xJ becomes 



Pa{x,(l)) = 



v2"+ix r 



Y exp(— M ) dw. 



(4.5) 



It can easily be checked that if \x\ > 1, the right hand side is bounded from above by Cx (^+2") and from below 
by cx~^^^'^"\ for some positive constants c, C. 



For the first derivative one can see that 



/ T T\ --y^au-\ ^exp(— a; ), 

{u^+x^)<^+2 ^/nx ^ ' 



which can be bounded using the same argument as before to obtain 

\p^^\x)\<C,x-^^+^-\ 
for some Ci and \x\ > 1. For the second derivative the procedure is the same. 

Lemma 4 Let r° be the diagonal matrix with elements 



Ther, 



' 47r7o l + a-isinh^/Ji 



dt. 



□ 

(4.6) 
(4.7) 



Proof. We analyze first the dependance on a for a fixed k. The functions {Lk{u)e "^^}fc>o form a orthonormal 
basis of L2(IR)- By using Cauchy-Schwarz followed by Lemma 5 from [S] we get 



l<J<- 



47r \Jo (1 + a-isinh2/3t)2 



1/2 



dt] < C (log af/'^, 



for some positive constant C. 
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Let now a be fixed and look at the asymptotic behavior of t^^ as A; ^ oo. We use the differential equation of 
the Laguerre polynomials, 12J 8.979: 



Thus 



Ln{x)^-iix-l)L'^X^)-xL';ix)). 
n 



-L„(^V2) = ti'„(^V2) 



L„(t72) = L;(t72) + iX(t72) 



which implies 



and 



1 1 d 

T^t -r- 

2 dt 



.-id 



2n 

Using integration by parts we obtain the formula 



dt 



dt^ 



1 

■dt= - 



Lk{t^/2)e-'"'^f{t)dt, 



1 + a-i sinh^/?t " k Jq 
where the function / is given by 

Pi{t) a-'^{P2{t)smh2l3t + P3{t)cosh2j3t) P3{t)a-'^ sinh^ 2 l3t 



1 + a-isinh^ f3t 



(l + a-isinh^ Pty 



(1 + a-i sinh^ /3i)3 ' 



(4.8) 
(4.9) 



with Pi(t) polynomials with degree at most four, whose coefficients do not depend on a. 

We split the integral into and and use the following bounds for the behavior of Laguerre polynomials in 
the two intervals (see [21] Theorem 8.9.12 and Theorem 7.6.4): 



max e-"/2|^„(x)| =0(n-i/4). 



and 



L„(x) = x-'/*0{n-'/*), 
uniformly on (0, 1]. Thus using Lemma 5 of [5], 



Air 



— / Lk(ty2)e- 



1 + sinh^ /3t 



■dt 



o((loga)4fc-5/4 



Lemma 5 For any (/3, L) there exists an a > such that Wa belongs to the class yV(/3, L). 



□ 



Proof. By using Minkowski inequality we get 

2 



\Jo 



r exp — r 



dw = 27r 

o l + z 







, -fa(z) 

ir- cxp — r 



2(1 -z) 



+ 2/3r dr 



1 - z 

1/2 



1 - z 



2 



dr < 



dz 



The interior integrals satisfies the bound 

r exp 



2(1 -z) 



-r'^rrh^^ + 2/3r ] dr < C(/3)(l - z). 
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for some positive constant C{/3). Thus 



Wciw) dw < C(/3) / a(l - z) 



1 X 2 



= C{fi) 



a + 1/2 



0, 



as a ^ 0. 



Lemma 6 For a < 1/2, the Fisher information of the family of densities TZW^ satisfies 

/(c)-C*(loga)X"'(l + o(l)) 

where Rq is defined in p.38p . 



□ 



(4.10) 



Proof. We sketch the proof foUowing the hne of 5_ and pointing out where the differences appear. After some 
transformations, the Fisher information of the family can be brought to the form 



/(c) 



1 



7^[lF-]([z,0]-w,, 



■ du. 



(4.11) 



By expanding 7?.[VF(f]([2;, 0] — u,(j)) ^ up to the second order and bounding the second derivative one can show 
that 



/ 



7^[M/-]([z,0]-u,, 



■ du 



Hl{u)du 



O 



u^Hl(u)dv^ =0{{\ogaf) 



as a — > oo. 



(4.12) 



RecaU that 'R[W^]{[z, 4>],cj))= 7^[^^=](0, cj)) does not depend on 0, cf. ([337| . Thus we can write 

1 



I{c) 



7^[M^c](O,0) 



Ht{u)du + 0{{\ogaY) 



From [5] we have 



/ H^(u)du = 1 



(4^/3) 



2(loga)3(l + o(l)) as a 



(4.13) 



From equations ()4.11|) - ()4.13p one obtains the desired resuh. 

The main difference in the proof appears when deriving the bound (j4.12p . To derive this, one needs a bound of 
the absolute value of the second derivative 



1 



dt'^n[w^]{t,i 



= o(i). 



as a ^ oo, 



t—\z.<h\—u 



which is uniform in u. We have 



where 



1 


< 

t—[z,(f>]—u 




+ 


P'ciuA? 


dt^n[w^'](t,(i)) 






Pc{u,(t))^ 



(4.14) 



Pc{u,cb) = n[wm^,4']-u,^)^n[Wc]{-u,cb)^n[Wc]{u,^)=Pa.{u,cb) + cHa{u), 

and the derivatives are with respect to the first argument. If w is in a compact interval, the two terms on the 
right hand side of (|4.14p are 0(1) as a ^ oo since in that case supj^j^,^^ \cHa{u)\ 0, and similarly for its 
derivatives. Outside the interval, we have that Ha and its derivatives are exponentially decreasing, and the 
dominating term is Pa- Indeed, according to Lemma[31 Pa is of order |it|~^^+^"\ and its first two derivatives 
are bounded from above by and respectively. 



□ 
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